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Abstract 

Within the real-time formalism (RTF) of thermal field theory, we apply the hard 
thermal loop (HTL) resummation technique to calculating effective two-loop thermody- 
namic potential in quark-gluon plasma (QGP) and its renormalization. The result with 
collective effects is obtained, which is valid for an arbitrary number of quark flavors 
with masses. 
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1 Introduction 

It is expected that at very high temperature and/or densities hadronic matter undergoes the 
deconfinement transition into a plasma phase composed of weakly interacting quarks and gluons 
(QGP). The basic thermodynamic characteristic of QGP system is its thermodynamic potential, 
and all other thermodynamic quantities of interest, such as the equation of state, can be derived 
from the potential by standard thermodynamic relations. In J.I. Kapusta's early works]!]], the 
thermodynamic potential of QCD at finite temperature was calculated by using bare propagators 
and vertices. However, without considering the influence from collective effects on propagation of 
partons in QGP, the result is not suitable for the QGP situation. Recently, P. Arnold and C. Zhai 
computed the free energy density for QCD at high temperature and zero chemical potential^, 
but for the case of massless quarks, and this technique what they developed seems difficult to deal 
self-consistently with various problems from the collective excitations in hot dense media. 

In quark-gluon plasma, deconfined partons are under hot dense condition and collective effects 
will affect their dynamic characteristic. For instance, interaction of gluons is no longer a long- 
range force but will change into a short-range ones. So when calculating thermodynamic potential 
'Supported by the National Science Foundation of China 



of the QGP, contribution from collective motion should be considered. On the other hand, R.D. 
Pisarski proposed the concept of hard thermal loop (HTL)|3| and from this a systematic method 
for the calculation of amplitudes in hot gauge theories is developed, which is termed the HTL 
resummation of E. Braaten -R.D. Pisarski[|]] and J. Frenkel and J.C. Taylor ||. Because the HTL 
reflects thermal fluctuations fairly well, in this scheme the physical effects induced by the collective 
motion in thermal medium can be discussed self-consistently. In this approach, it is natural to 
distinguish between " hard" momenta, ~ T, and the "soft" momenta, ~ gT. As a consequence, 
for a consistent perturbative calculation, soft momenta lines and vertices should be replaced by 
their effective counterparts in which the HTL is introduced. It is well-known that thermodynamic 
potential is the generating functional of vertex function at zero momenta. Obviously, its external 
lines is all soft, hence when performing the calculations of the potential the HTL should be taken 
into account. 

In this paper, we will evaluate the effective two-loop thermodynamic potential with Nf massive 
quarks based on the Braaten-Pisarski technique, and we give a effective method to deal with the 
renormalization of this effective perturbation theory. We work in the real-time formalism. In order 
to handle the ultraviolet singularities which occur in our procedure, dimensional regularization 
(DR) is used. The D dimensional vector, P^ = (p°,p), is contracted with Minkowski metric, 
P 2 = p 2 — p 2 and p =| p |. For a internal momentum P, we introduce the following notation: 

d D P 
(2tt) d 

The plan for the rest of the paper is as follows. Section 2 gives a description about HTL gluon 
self-energy and its effective propagator. Sec. 3 computes the effective thermodynamic potential with 



f (dP)s f 



fermions in RTF and the renormalization. Sec. 4 gives a summary and conclusions. 

2 HTL gluon self-energy with massive fermions and gluon effec- 
tive Propagator 

Since the momenta of partons in QGP can be soft in our considerations, we will need to 
use effective propagators including the HTL corrections. In covariant gauge the effective gluon 
propagator is decomposed into its transverse, longitudinal and gauge components Q: 

- & U {K) = A^(K)A T {K) + B^(K)A L (K) + ^ U {K)A G (K), (1) 

where £ is the gauge parameter, A^ v ', B^ v and D^ v are the transverse, longitudinal and gauge 
projectors, and 



where Ht,l{K) stand for two independent components, transverse and longitudinal terms, of the 
gluon self-energy Ii^ v {K). 

In the rest-frame of the QGP, using the projectors 

A 00 (K) = A i0 (K) = A 0i (K) = 0, A ij (K) = g ij + — 

k 

B^{K) = g^ -A» V {K)-D» V {K), D^(K) = ^-^-, (3) 

and the identity 

W(K) = -A^ U (K)U T (K) - B^ U (K)U L (K), (4) 

we have 

uiuj 

(D-2)U T (K) = -(gV + —)Il tJ (K), 

u L (K) = -(i-M)n 00 (iO, (5) 

where D = 4 — 2e denotes the number of space-time dimensions. 




Figure 1: One-loop gluon self-energy. Curly lines, solid lines, and dashed lines represent the 
propagators of gluons, quarks, and ghosts, respectively. 

The gluon self-energy at one- loop level is given in Fig.l. The explicit form of the bosonic 
contributions, graph (a), (b),(d) of Fig.l, Ti T L (ko,k) are[f?J 

n?(k ,k) = Sm^Kl-^-^^ln^)], (6) 

where x = -£ and rn} g = Ng 2 T 2 /9 describes the screening due to thermal pure gluons. 

For the fermionic case, some results obtained come from the massless theory only[^, |7j. Here, 
we derive the 11^ L (ko,k) with fermions masses in the HTL approximation. We note that the 
physical component of the fermion propagator at finite temperature has the form 

S F (P) = (P + m)[ p2 - 2 - 2ir5(P 2 - m 2 )n F (P)} 
P L — m z 

= (P + m)(A° F (P) + A F (P)), (7) 



3 



where n F (P) represents the fermion distribution function. The temperature dependent contribution 
of the fermion loop, graph (c) of Fig.l, to the self-energy is given by 

- m$ F (K) = N F T F g 2 J (dP)T r [ 7M (P+ JQ 7 „ P ][A° F (P + K)A F + A F (P + K)A° F (P)} 
= -2mDT F N F g 2 [ ( dP )^Ef--^l[(P . K){P,K U + P U K,) - K\P,P V ) - g, v {P • K)\ (8) 

j {p-Ky-(—y 

where the group theory factor T F 5 ab = T v (T a T b ), N F is a number of flavors and the terms with 
two statistical factors has been dropped by 5-function constraints. From Eq.(5) and (8) one obtains 



D-2 J K ' ( P . K )2 _ (Kly a ; fc2 w 

Ul{K) = 2nDT F N F g> J (dP) "*™^ ^ [{P • + ^ - 2poW • if)]- P) 

Now we apply the HTL approximation. Let the external momentum K be soft, ~ gT, and the 
internal ones P be hard, ~ T, we can approximate 



2 - 



(p . #)2 _ (2|i)2 (P-if^ 4(P-if) 

= (7r^(i + 0(5 2 ))- (io) 

Inserting Eq.(lO) into Eq.(9) and integrating over po and 9 (p • k = pk cos 9), we get the 
fermionic result with masses: 

U T (K) = [ d ppv[x* + (1 - x 2 )^- ln(^)]n F (E), 

7T J ZV X — V 

U F L (K) = [dppv(l-x 2 )[l + 4^--H—)]nF(E), (11) 

7T J X z — V z V X — V 

where E = \/m 2 + p 2 , and v = %- is the velocity of the quarks. Using Eq.(6) and Eq.(ll), one can 
obtain the HTL gluon self-energy with N F massive flavors, 

U T , L (K) = Ii TL (K) + Iil L (K). (12) 

Because the poles of G^ U {K) correspond to eigenmodes of the QGP system, Eq.(l) and (12) 
gives us the desired spectrum of gluon collective excitations. For the sake of describing the screening 
phenomenon of thermal gluons, the screening mass is introduced and defined as|| 

Ml L = lim II T , L (0,fc) (13) 

k 2 ->-M%, L 

According to this definition, we have 

U T (k = 0,k) = 0, 
U L (k = 0,k) = M 2 g 

= \Ng^ + ^fjdp{E + ^)n F {E). (14) 



It implies static screening of the chrome-electric field but not the chromomagnetic field. In order to 
avoid being in infrared troubles, a magnetic mass for the transverse gluon is introduced by assuming 
that it modifies the static behavior in such a way thatQ 

n T (o) = Ml, (15) 

where M m ~ g 2 T. 

Substituting Eq.(12) and (13) into Eq.(2), at finite temperature, one can obtain the components 
of the HTL gluon effective propagator: 

A t,l( K ) = K 2_ l M 2 +^{K 2 -M 2 m Jn B {K) 

= A° TjL (K)+A^ L (K), 
A G {K) = + 27r5(K 2 )n B (K) 

= A%(K) + A G (K), (16) 

where n B {K) represents the boson distribution function. 

3 Effective thermodynamic potential in RTF and its renormaliza- 
tion 

In this section, we determine the soft exchange corrections to the thermodynamic potential 
in QGP , the relevant diagram to be analyzed is depicted in Fig. 2, in which the exchanged gluon 
has momentum components much smaller than T while the fermions and other gluons have hard 
momenta. According to the idea of the HTL resummation, for the soft exchanged gluon effective 
propagator G^ U (K) must be used, but the bare ones are needed sufficiently for the others, and in 
our situation no effective vertices is needed because the vertices are connected, at least, to one hard 
line. 




Figure 2: Effective two-loop thermodynamic potential. The solid blobs indicate the HTL effective 
propagator. 
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Using the standard Feynman rules of QCD to evaluate the diagram of Fig. 2(a), one obtains 

I a = -C A g 2 ^ j \dK)(dP)[-(2K -P) x9llu + (K -2P) v g Xll + {P + K) lx g vX ][-g rv + {l-Z)D TV (K)] 
xA G (K)[-(2K - P) p g Ta + (K - 2P) T g ap + (P + K) agpT ][-g Xp + (1 - £)D A '(P)]A G (P) 
x [A^(K - P)A T (K -P) + B pa [K - P)A L (K - P) + £D^(K - P)A G (K - P)\ 
= lP+li L) +lP, (17) 

where the group theory factor C A 5 ab = f acd f bcd , n is a mass scale in the DR and 

If) = CWV e / (dK)(dP)[(4D - 6)P^P U A^(K - P) + (D - 2){(2K - P) 2 + (2P - K) 2 )} 
x A G {K)A G (P)A T (K - P), (18) 

if) = CWV e J(dK)(dP)[(AD-6)P p P l/ B^(K-P) + ((2K-P) 2 + (2P-K) 2 )} 

x A G (K)A G {P)A L (K - P), (19) 

J«) = C A g 2 ^ J(dK)(dP)[(AD-Q)P p P u (g^-D^(K-P))-(D-l)(K + P) 2 ] 

x A G (^)A G (P)A G (K-P) (20) 

are the transverse, longitudinal and gauge terms, respectively. 

Substituting Eq.(3), (14) into Eq.(15), one can evaluate these integrals by the dimensional 
regularization method at finite temperature. Since the real part of thermodynamic potential is 
concerned, those terms with three or no statistical factors will be dropped [p7c|] . In Feynman gauge, 
including one distribution function term, Iip ta , and two distribution functions term, l2/3, a , we arrive 
at 

hp,a = -iC A g 2 ^ej^^{(f + lle)(M 2 m )^J(dK)A^K) + (f + 8e)(M 2 )^ j (dK)A? L (K) 

+5[2(M 2 J 1 -* + (M^) 1 -] / (dP)A G (P) + 10e J (dP)A G (P)(p 2 )^[F(0,M 2 ,p 2 ) - F(0, M^p 2 )]}, (21) 

hp,a = C A g 2 |(dK)(dP){8M^(2A G (^)A G (P)A^if-P) + A°(K-P)A G (P)A G (K)) + iA G (P) 

x [9(A G (K) + 2A? L (K - P)) + 4(A?(X - P) - A{{K - P)) + 7{A G {K - P) - a{{K - P))] 
+^M 9 2 (2A G (K)A G (P)A^(^ - P) + A° L (K - P)A G (P)A G (K))}, (22) 

where function F(m 2 , M 2 ml p 2 ) is defined as 



F(m 2 , Ml m ,p 2 ) = dy {l v) * {[M 2 ^ - m 2 + ^ \M* >m - m 2 ) 2 - \yp 2 M* >m ] 

P J o y 



M 2 . m -m 2 - 2yp 2 + J (M| - m 2 ) 2 - 4yp 2 M 2 m 
x In — — : : ^=L = — + [M 2 m - m 



M 2 m - m 2 + J(M 2 m - m 2 ) 2 - 4yp 2 M 2 



- J M 2 - m 2 2 - 4yp 2 M 2 x In v = }, (23) 

where y is a Feynman parameter. 

It is easy to find that l2/3,a is finite. However contains, in addition to finite parts, ultravio- 
let divergences which take the form of multiplying statistic factor. Further more, these overlapping 
divergences are thermal mass dependent, which differs from naive perturbation theory. This sit- 
uation needs a new treatment of the renormalization, and hence we adopt the counterterms in 
Fig.3(a), (b). Where the parentheses indicate that enclosed subgraph is evaluated in the vacuum 



(a) (b) (c) (d) 

Figure 3: Renormalization counterterms including three-point gluon vertex and four-point gluon 
vertex 

limit, and the particle lines entering and leaving the subgraph are put on shell. 

Because of using HTL effective perturbation theory, when performing the procedure of sub- 
tracting divergence, in addition to zero-temperature subgraph Fig.3(a), one must consider the 
contribution from a thermal subgraph, like Fig. 3(b). And when calculating the vacuum limit of 
a subgraph, one should save the corrections in propagator and vertices from the collective effects 
induced by HTL, and discard the terms with thermal distribution function. 

Considering the Fig. 3 (a), one can get the pole term of the subgraph easily 

K = iC *9 2 (i^t-y^ + y^Vl- (24) 

For the Fig.3(b), we may multiply by g^ u to simplify the algebra, the relevant part of the subgraph 
is 

g^lgg = C A g 2 ^ J (dQ){(4D - 6)P x P p A xp (Q)[A^(Q) - A° L (Q)} 

+ [(P + 2Q) 2 + (Q- Pf] [(D - 2)A° T (Q) + A° L (Q)]}A° G (P + Q). (25) 

Inserting projector A Xp (Q) into the above equation and with a bit of operation, we find that the 
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first term is finite, and applying mass-shell condition to momentum P, we have the pole term 

g^vl = i5C A g 2 j^- e [2Ml + M 2 }. (26) 

When these pole terms are inserted back into the diagram Fig.3, one obtains the counterterms of 
the overlapping divergences of the Fig. 2 (a) 

lT V) = J(dPKG^(P) 

= -^A9 2 ^J(dP)[fMiA^P) + fM 2 A^P)}, (27) 

iT v) = - j \dP)vi g9 ^ G {p) 

= ^C A g 2 j^[2M^ + M 2 } J (dP)Ag(P). (28) 
So, after considering the renormalization, total ultraviolet divergences of the Fig.2(a) is 

j{div) r j{div) M{div) T P{div)i _ n fr>Q \ 

finite part is 

A fin) =h P ,a + A f $, (30) 

The steps to evaluate the Fig. 2(b) is similar to the above mentioned, we obtain 

h = -2{D-l)iC A g 2 ^ J(dQ)(dP)[(D-2)A T (P) + A L (P)+A G (P)]A G (Q) 

= h/3,b + h/3,b, (31) 



where 



l w>b = -2iC A g 2 ^ ^±^ [(6 - lOeJCO 1 " 6 + (3 - 2e)(M 2 )^] j (dQ)A G (Q), (32) 
■ 2/3)6 = -2i(D - l)C A g 2 ^ J(dQ)(dP)[(D - 2)A?(P) + Aj(P) + A G (P)]A G (Q). (33) 



The relevant counterterms are shown in Fig. 3(c) and (d). The graph (c) of Fig.3 gives zero, 
since its zero-temperature subgraph vanishes in dimensional regularization. For Fig. 3(d), after 
straightforward and tedious calculating, thermal subgraph yields pole term 

Vf g = -2iC A g 2 j^\ gixv {2M 2 m + \m 2 ) - g»g»g l3 2 -{M 2 m - M 2 )], (34) 

therefore, the subtractive terms are 

llf v) = 0, (35) 



8 



J (dP)vf g g^A G (P) 



= 2iCV^^[(6-^)M^ + (3-^ e )M 9 2 ], 
then, total overlapping divergences and finite part from Fig.2(b) read 



(36) 



(div) _ 1 r Jdiv) T 0(div) T /3(div) 



0, 



(37) 

4 /m) = /2A6 + 4£ ) -4 (/m) - (38) 
As for a ghost loop in Fig.2(c), which its counterterms are shown in Fig.4(a), (b), one can get 

I c = C A g 2 ^ J (dK)(dP)K^P v A g (K)A g (P) 

x [A^A T {K -P) + B^(K - P)A L (K - P) + A^(K - P)A G (K - P)] 

= h/3,c + h/3,c, (39) 

where A g (K) is a ghost propagator, and 

hp,c = -iC A g 2 ^^^{\{l-e){M 2 m f^ J{dK)A^{K) 



(4vr) 2 

+l(l + 8e)(M 2 ) 1 ^ f{dK)A{{K) 



+ J (dP)A§(P)(p 2 ) 1 ^[F(0,M^,p 2 )-F(0,M 2 ,p 2 )}}, (40) 

C A g 2 ^J (dK)(dP){ l -[A§(P)(A G (K - P) - A^K - P)) + 2A^P)A^ L (K - P) 
M 2 

-Af (If)Aj(P)] - -f[2A%K)Al{P)A{{K - P) + A° L (K - P)A^(K)A^P)} 

+ b 2 " (P (k (k p) P 2 ))2 ][(Ai(if - P) - A° T (K - P))tf g (K)*P g {P) 

+2(A ( l(K - P) - A?(K - P))A° g (K)AP g (P)}}. (41) 




(a) 



(b) 





(c) 



(d) 



Figure 4: Renormalization counterterms including ghost-gluon vertex and fermion-gluon vertex 
The pole terms and counterterm contributions from Fig.4(a), (b), respectively, are 

Kg = tC ^ 2 J^T e [\K,K v + ±K 2 g, u ], (42) 



T 0(div) 
L gG 



Consequently, we obtain 



V & = %CAg2 2(^fe K ^ (43) 
J (dK)V g ° G G^(K) 

l ° A92 Wfe I ^^ M M^) + \mIK(K)], (44) 

= J (dK)V&*&{K) = 0. (45) 



T (div) _ 1 Adiv) 1 Mdiv) T P(div) _ „ , 

I { J m) = Iw, c + I% m J. (47) 



For a soft momentum (i^ — P) ~ gT, we have 



"»^-^ = 5TC^rr7 "iA^T- (48) 

Substituting Eq.(46) into Eq.(28),(36) and (45) with a little practice, the contributions, up to <? 3 , 
to the thermodynamic potential from bosonic graphs, (a),(b) and (c) of Fig.2, together may be 
written 

I B {g,M g {m),T) = ±[#*0 + + 

= iC A g 2 — [— + ^-(9 + —Ml aTctan EL)] +0 (g 4 ). (49) 
y 24 L 24 2vr 2V 2 gT M g vy ' y ' 

The bosonic sector has been evaluated from the above. Now, we shall incorporate fermions 
into the theory. The diagrams to calculate are as shown in Fig.2(d) and Fig.4(c), (d). The loop 
integrals with fermions can be handled in similar manner as bosonic loops. For graph (d) of Fig.2, 
we write 

I d = -T F N F g 2 ^ J (dK)(dP)T r [^S F (K) lu S F (P)]G^(K - P) 

= hp,d + l2/3,d, (50) 

here, what remains is the part with one or two statistic factors. Performing the trace, applying the 
Eq.(3) and the dimensional regularization, after a lengthy calculation one can get 

hp,d = -iT F9 2 ^f(4 - 2 e )^|(^)1m 2 / (dP)Aj(P)[3 + ,(14 - M In ^ 
M , M Q ~ \/ M n ~ 4m 2 m 2 M 2 M 2 M 2 



m 2 V 9 mm /m2 _ j™2 m 2 m 2 m 2 m 2 



M g + ^M 2 - 4m 2 
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5 m 2 1 / 2m 2 — M 2 + M g J M| — 4m 2 

x [1 + e(? " Stt^t + -z^rJMt-lm 2 In \ g )] 

L l 3 M| 2M g V ° 2m 2_ M 2_ M ^ M 2_ 4m 2^ 

+ ^|(dP)A^(P)/[F(m 2 ,M 2 ,p 2 )-F(m 2 ,M^,/)]}, (51) 
/ 2/M = T F g 2 ^f{D) J (dK)(dP){i(D - 2)[A F (K)A F (P) - (A?(tf - P) + A^(K - P)K(P)] 

+ (2m 2 + ^M 2 )(A° L (^ - P)A£(K)A£(P) + 2A F (K)A^(K - P))A?(P)) 
+^M 2 (A° (K - P)A£(^)a£(P) + 2A F {K)A? T {K - P)Aj(P)) 
+ 2 ^ 2 " (P (k (k p) P 2 ))2 ][(A° L (^ - P) - A° (K - P))A F {K)A F {P) 

+ 2(A^ - P) - A^(K - P))A° F (K)A^(P)]}. (52) 
The pole terms and counterterm contributions from Fig. 4(c), (d), respectively, are 

V° G = -iT F g 2 j^-(-^K,K u + ^V). ( 53 ) 

Vj* G = -*(T a T a ) Q/35 2 ^-^(P - 4m), (54) 
= j (dK)V^ G G%(K) 

= - iT F9 2 j^ e \J{dK)[2Ml^ T {K) + M 2 g A{{K)l (55) 

= - 1 (dP)T r [F/ G (P)(P + m)]Aj(P) 
= _*3m 2 T F5 2 ^^/(4-2 e ) j(dP)A£(P). (56) 

As a result total ultraviolet overlapping divergences of the Fig. 2(d), with considering the renor- 
malization, is 

T (div) _ 1 Jdiv) 1 T 0(div) T fS(div) _ „ /c^A 

J d - 2 i/M ~ 2 FG ~ FG ' ^ ' 

To compute the contribution from finite parts, we take the limit e — > 0. Integrating over po and fco 
and changing the variable of integration by letting (k — p) — > k, we get 

n F (E p )n F (E k ) f f n F (E p ) n F (E p+k ) 



Mdiv) 
FG 



T (3(div) 
FG 



t o rr 2 /0 /* /. ^r^ n F{Ep)n F {E k ) f f n F (E p ) i 

hf3, d = 2iT F g {2 (dp)(dk)—- +/ (dp) (dk){—= — [- 

Jp,k(hard) &p &k Jp(hard) Jk(soft) &p 



E 



p+k 
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x [(2m' + M + 2(p 2 - ^-))( t= i t-t -To ro + 



(P • k) 2 

9 k 2 "\E p + £ p+fe ) 2 - fc2 _ M 2 < {Ep _ Ep+k) 2 _ k 2_ M 2 g 

+ (M 2 - 2(, 2 - ^))( __^_ + __^__ )]] - [^1(58) 



x(2m 2 + M 2 + 2(p 2 - ^))(-___^. + 



k 2 ' ( jBp + £ ;f))2_( k + p)2_ m 2 ( £p _^))2_ (k + p) 2_ m 2 

+ ""7* ' (M^^-^-fHX m + w )]}}, 



«B(4 T) ),^2 ,2 (p-k) 1 



where 



When momentum fc is soft, we can approximate 



Ek+p + E p m 2E P H — — , Ek +p — E p — — , 

= e/?£ J p - x ~ nF(^p). (59) 

Substituting Eq.(57), (46) into Eq.(56), after a tediously long calculation we obtain, up to g 3 , 
contributions from the loop with massive fermions 

lr. .(/in). 



I F (g,m,M g (m),T) = -[h M + I?& 



where 



iT F N F {2g 2 I 2 (m,p) - —— 5 3 [2m 2 (l - ^axctan f-)/ 2 (m, /3) 

{2tt) z gl M g 

-< x + \€; + %^§)hw + |$ + \§)^m (60) 

+^?V 2 ((f^) 2 - \{\ - ^arctan||))/ 5 (m,/3) +/ 1 (m,/3)]} + 0( 5 4 ), 



lim ' W " (2vr)2i *v^T? e ^ +1 ' 

2 /-oo ^2 y 

/2(ra,/3) = Wio dp m2+p2 (^^^ + 1)2' 

h{m ' P) = j2^J dP m 2 +p 2 (e /V^* + i)2' 

2 Z" 00 p 3 1 

h(m,/3) = / dp = r^-^ • 

(27r) z Jo V m +P (e^v™ +p + l) 2 

2 Z" 00 p 2 1 



%^) = 7^I fr , / ... -j == • (61) 

(z7r) jo (m 2 +p2) 5 e /3V m2 +P 2 + 1 
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Taking some proper approximations, such as a high-temperature expansion, these integrals in 
Eq.(59) can be evaluated. 

Combining Eq.(58) and (47), finally we get the effective two-loop thermodynamic potential in 
QGP 

tt exch = -i[I F (g,m,M g {m),T) + I B {g, M g (m),T)], (62) 

in which the effects both of finite fermion masses and of collective motion induced by the HTL are 
involved. Using this expression for the thermodynamic potential, we can consider the limiting case 
of our interest. 

In the ultra-relativistic (m = 0) limit, Eq.(60) yields 

= ^(C, + 2T FWF)+ /(^ {[9 + f^,A f + f)a I c tav /Z|Z 1 £| 

r 1 1 /, 2 2 / 6 FT~ T iV^7. / 6 2.. 

+ [ 2 + ^ (1 + 2FT^ " 3 V ^ + V 3 (iV + — )(arCtan V WTl^ - 3 )} 

x(y + ^C(3)-2)]T F iV F } + 0(< 7 4 ). (63) 

This result agrees with the one found in the imaginary time formalism (ITF) 0, which is to be 
expected because of our considering only HTL in the static limit. 

4 Summary and conclusions 

In the present paper, we have derived the HTL gluon self-energy with fermion masses and its 
effective propagator by use of the RTF. we have calculated explicitly effective two-loop thermody- 
namic potential in QGP, and have made a detailed presentation of its renormalization. we obtain 
the result including collective effects due to the HTL, which is suitable for an arbitrary number of 
quark flavors with masses. 

In particular, we have found that the effective two-loop thermodynamic potential in QGP is 
composed of two parts. One is the contribution, g 2 terms, from interactions between the hard bare 
partons, which is just the result in naive perturbation theories. The other part is the higher order 
terms, which is owing to the interactions between soft gluon and hard partons in hot dense media, 
and chromomagnetic mass effects is absent from the g s level. 

To be sure, after introducing the THL, one can still apply the DR to handle self-consistently the 
regularization and renormalization of the thermodynamic potential, but a point which is especially 
to be noted is that, within this effective perturbation theory scheme, we should introduce the 
thermal counterterms and employed the RTF, which makes our calculations simple and clear. 
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Finally, it will be straightforward to evaluate the thermodynamic potential to higher order 
making use of the methods and formulas developed in this paper. 
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